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Abstract 

It is known that a 2 x 2 quaternionic matrix has one, two or an infinite 
number of ieft eigenvalues, but the available algebraic proofs are difficult 
to generalize to higher orders. In this paper a different point of view 
is adopted by computing the topological degree of a characteristic map 
associated to the matrix and discussing the rank of the differential. The 
same techniques are extended to 3 x 3 matrices, which are still lacking a 
complete classification. 
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1 INTRODUCTION 

In 1985, Wood proved that any n x n quaternionic matrix A has at least 
one left eigenvalue, that is a quaternion A £ H such that the matrix A — A Id is 
singular. However, even for matrices of small size the left spectrum is not fully 
understood yet (see Zhang's papers [Ml El] for a survey). For instance, it was 
only in 2001 when Huang and So |13j proved that a 2 x 2 matrix may have one, 
two or an infinite number of left eigenvalues; a different proof was presented by 
the authors in ^T\. While Wood used topological techniques, namely homotopy 
groups, the two latter papers are of algebraic nature and seemingly difficult to 
generalize for n > 2. 

In this article we try a different approach. The basic ideas will be those of 
characteristic map, linearization and topological degree. In the simplest case 
n — 2 we associate to each matrix A a polynomial iia whose roots are the 
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left eigenvalues; computing the rank of its differential allows us to classify the 
different types of spectra. 

In the second part of the paper we extend those techniques to 3 x 3 matrices. 
This time the characteristic map may not be a polynomial but a rational map, 
usually with a point of discontinuity. This seems to require the use of a local 
version of degree, although a closer look allows us to reduce the problem to 
the global theory. In particular, this gives a new proof of the existence of left 
eigenvalues. 

Unlike the 2x2 case, where the linear equations that appear correspond to 
the well known Sylvester equation, the 3x3 situation is much more complex. 
Then for n = 3 a complete classification of spectra is still unknown. Nevertheless 
our method allows to deal with specific examples and opens the way for a better 
understanding of the general case. 

The paper is organized as follows. In Section [5] we recall some topological 
and algebraic preliminaries. Although our ideas are closely related to the theory 
of quasideterminants |10j . we have preferred a development based on Study's 
determinant which parallels the commutative setting. Section [H] is devoted to a 
notion of characteristic map for the left eigenvalues of a quaternionic matrix A, 
that is a map /i^ : HI — >■ H such that /za(A) = if and only if the matrix A — Xld 
is not invertible. In Section|3]we give a complete classification of the left spectra 
of 2 x 2 matrices, depending on the rank of the characteristic map. In Section [S] 
we prove that any 3x3 matrix A has a characteristic map fj,A which is either a 
polynomial (when some of the entries outside the diagonal is null) or a rational 
function with a distinguished point tta called its pole. When fiA is continuous 
it has topological degree 3. However, tta may be a point of discontinuity; in 
this case the matrix B = A — ttaU turns out to be invertible and we prove 
that B and B~^ have diffeomorphic spectra and that B^^ admits a polynomial 
characteristic map. The last Section offers several illustrative examples. 

2 PRELIMINARIES 

2.1 Topological degree 

The topological degree (or Brouwer degree) of a map can be defined by tech- 
niques either from algebraic topology [5] or from functional analysis (TJ [H [23] . 
We want to apply the following global result (cf. [19l p. 101]): 

Theorem 1. Let M be a connected closed oriented manifold, let fi: M ^ M 
be a differentiable map of degree k. Let m ^ M be a regular value such that the 
differential preserves the orientation for any A in the fiber ^~^(rn). Then 
fj,~^{m) is a finite set with k points. 

Sometimes one has to deal with a local notion of degree. The main result is 
the following one [I] p. 38]. 

Theorem 2. Let Q. be a bounded open set in R". Let /i: — > M" be a map 
which is continuous on the closure and differentiable on Vt. Suppose that is 
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a regular value of ji and that ^ /i(9f7). Then 

A6ai-i(0) 

where we denote by Jf^ the Jacobian of fj,. 

A well known consequence is that if dcg(/i, £7,0) ^ then the equation 
^(A) = has at least one solution in fl. In fact, for maps of the sphere into 
itself, the existence of solutions only depends on the continuity, by the following 
result [201 Ch. VIII, Ex. 2.5, p. 191]: 

Proposition 3. Let fi: S"^ ^ S'^ be a continuous map whose degree is nonzero. 
Then fi is surjective. 

2.2 Linearization 

We consider the space of quaternions H as a difFerentiable manifold diffeomor- 
phic to R"^. Then the differential : H ^ H at the point A G H of the 
difFerentiable map /i : H — H can be computed by means of the formula 

M*a(X) = ^ ^(A + tX) - hm ) (/i(A + tX) - /i(A)). 

at \t=o t^o t 

Lemma 4. 1. Let /, g: H — )■ H be two differentiable maps. Then the differ- 
ential of the product is given by 

if ■ g%xiX) = • .9(A) + /(A) • g,x{X); 

2. assume that /(A) 7^ for all A G H. Let 1//: H ^ H be the map given by 
(l//)(A)=/(A)-i. Then 

{l/fU{X) = -f{\)-'f^,{X)f{X)-\ 

2.3 Sylvester equation 

Let P, Q, i? G H be three quaternions. We are interested (see formula ([6])) in the 
rank of R-linear maps E : H ^ H of the form E(X) = PX + XQ. The equation 
E(X) = R has been widely studied, sometimes under the name of Sylvester 
equation |17) . 

Lemma 5. 1. Let P = t + xi + + zk. Then the matrix associated to the 
R-linear map X ^ PX with respect to the basis {l,i,j,k} is L{P) = 
5ft(P) Id + A{P), where 3?(P) is the real part of P and 



A{P) = 
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2. analogously, if Q = s + ui + wj + wk then the matrix associated to the right 
translation X ^ XQ is R{Q) — Id + B{Q), where 



B{Q) = 



Next Theorem is a reformulation of the resuhs by Janovska and Opfer in 
[T5] . see also [TT] . 
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Theorem 6. 



i. The rank of T, is even, namely 0, 2 or 4; 



2. rankE < 4 i/ and onZ?/ if P and — Q are similar quaternions, i.e. 
have the same norm and the same real part; 

3. ranks — if and only if P is a real number and Q = —P. 
Proof. The matrix associated to S is 



they 



J = 



t + s 
X + u 
y + v 
z + w 



—X — u —y — V 
-z + w 
t + s 



t + s 
z — w 
—y + v 



X — u 



—z — w 
y-v 

—X + u 
t + s 



Since 



det J ={t + sf + 2{t + sf{x'^ + y"^ + z"^ + u"^ + v'^ + w^)+ 



(1) 



(.2 



y 



the matrix J has rank 4 excepting when t+s = and x +y +z = u +v +w . 
In this case J is skew-symmetric, hence its rank is even. If rank S = it follows 
that X — y = z = and u = v = w — 0. □ 

2.4 Resolution of arbitrary linear equations 

More generally, let us consider a linear equation of the form 

PiXQi + --- + PnXQn = R, with P„Q,,i? e H. (2) 

For each bilateral term X M- PXQ, the matrices L{P) and R{Q) commute, 
because P{XQ) = {PX)Q. Then A{P) and B{Q) commute too. This implies 
that the quaternionic linear equation ([2]) is equivalent to the real linear system 
MX = R, where X,ReR'^ and M is the 4 x 4 real matrix 

M = J2 L{P,)R{Q^) = <^^b^ Id + 51 («^^» + ^^^0 + J2 

i i i i 

with ai = 5R(F,), b, = ^(Qi), A, = A{Q,) and B, = B{Q,). Contrary to the 
case n — 2, when n > 3 the rank of M may be odd. 

Example 7. The rank of the matrix associated to the bilateral linear equation 
kX + X{2 - i) - 2jXj = is 3. 
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2.5 Determinants 

It is possible to generalize to the quaternions the norm | det | (that is, with 
real values) of the complex determinant. Papers [H |3] are surveys of the gen- 
eral theory of quaternionic determinants. For the relationship between Study's 
determinant and quasideterminants see |101 pp. 76-85]. 

Definition 8. Let the quaternionic matrix A G A^„xn(EI) be decomposed as 
A = X + jY with X,Y e A4„xn(C). We shah call Study's determinant of A 
the non-negative real number 



Sdet(A) := (dctc(A))i/2^ 



where c{A) is the complex matrix c{A) 



X -Y 
Y X 



e M2nx2n{C). 



Remark 9. Up to the exponent, this is the same determinant which appears 
in Theorem 8.1 of [26], as well as others considered in [2j. We have normalized 
the exponent to 1/2 in order to ensure that Sdet(_D) = \qi . . .qnl for diagonal 
matrices D = diag(gi, . . . g„). 

Proposition 10 ([3j). Sdet is the only functional that verifies the properties: 

1. Sdet(AB) = Sdet(^) • Sdet(B); 

2. if A is a complex matrix then Sdet(A) = | det(A)|. 

The following immediate consequences are very useful for computations. 
Corollary 11. 1. Sdet (A) > if and only if the matrix A is invertible; 

2. let A and B PAP^^ be similar matrices, then Sdet(^) = Sdet(_B); 

3. Sdet (A) does not change when a (right) multiple of one column is added 
to another column; 

4- Sdet (A) does not change when a (left) multiple of one row is added to 
another row; 

5. Sdet (A) does not change when two columns (or two rows) are permuted. 

We shall need the following result too (we have not found it explicitly in the 
literature): 

Proposition 12. For any matrix with two boxes A, B of order m and n respec- 

^ ^ = Sdet (^) • Sdet (B). 



tively, it holds that Sdet 
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2.6 Jacobi identity 

Let C be a complex n x n matrix. Let / — {ii, . . . ,ip} and J = {ji, . . . ,jp} 
be two subsets of {1, . . . ,n} with the same size p. Let us denote by Cjj the 
submatrix formed by the rows with index in / and the columns with index in 
J. On the other side, let us denote by C^ "^ the complementary matrix obtained 
by suppressing the rows in I and the columns in J. 

The following Jacobi identity is attributed to Kronecker in [2 2 j . 

Lemma 13. Assume that the complex matrix C is invertible. Then 

det {C-^)ij = (-1)^+-^ det C^^V dot C, 

where I + J means ii + ■ ■ ■ + ip + ji + ■ ■ ■ + jp- 

A generalization to quasideterminants appears in [10; Theorem 1.5.4, p. 74], 
see also Section IX2l We shall use Study's determinant to establish an analogous 
result in the quaternionic setting. 

Proposition 14. Let A be an invertible quaternionic matrix. Then 

Sdet {A-^)i,j = SdetA-^'VSdet A 

Proof. If / = {ii, . . . ,ip} we denote /' = I+n = {ii+n, . . . ,ip + n}; analogously 
J' = J + n. The result follows from Lemma [T51 because 

c{iA-')j,j)=c{A~')jui',j.j' = (c(A)-^),,,, ,,,, 

and 

c{A-^'')=ciAy'^''''"^''. 

□ 



3 CHARACTERISTIC EQUATION 

The problem we are proposing here is to find a characteristic map for the left 
eigenvalues of a given matrix A, that is to find a map fXA : HI — > 11 such that 
/Xyi(A) = if and only if A is a left eigenvalue of A. Notice that the function 
Sdet(yl — Aid) is real- valued, so it is not of interest from the point of view of 
the topological degree, nor it is solvable in any obvious way. 

3.1 Left eigenvalues 

Let A he Si matrix with quaternionic coefficients. 

Definition 15. The quaternion A £ H is a left eigenvalue of A if the matrix 
A — Aid is not invertible, or equivalently Sdet(74 — Aid) = 0. 

Let <7i{A) be the left spectrum, i.e. the set of left eigenvalues, of the matrix 

A. 
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Proposition 16. The spectrum (Ti{A) is compact. 



Proof. The spectrum is a closed set because it is given by the equation Sdet(j4 — 
Aid) = 0. It is bounded because A e cr;(^) if and only if there exists v G H", 
V ^ 0, such that Av = Xv; then 

I ^ I l^^^l \AW\ , , , 

|A| = ^ < sup ^ = A . 

M \w\=i \w\ 

□ 

Proposition 17. Let B be an invertible matrix. Then A G (ti{B) if and only if 

Proof If Bv = Xv then B-'^{Xv) = B-^Bv ^v = X-^{Xv). □ 
3.2 Background 

When the matrix A is hermitian, all left eigenvalues are real numbers so they 
coincide with the right eigenvalues [7]. Moreover it is possible to define a true 
determinant for hermitian matrices [181 [26] , which allows to construct a char- 
acteristic polynomial p{t) = det{A — t Id) with real variable. 

For the general case, Sdet(A) equals, up to an exponent, the determinant of 
AA* . On the other hand, Zhang |26| pointed out that if the quaternionic matrix 
is decomposed as A ^ X + jF, with X,Y & Mnxn{C), then its left eigenvalues 
X = X + jy, with x, y e C, are the roots of the function a: C x C R given by 



a{x, y) = det 



X~xld -Y + ylA 
Y~yld X-xld 



(3) 



This is equivalent to the equation Sdet(A — Aid) = 0. 

Another approach is due to Gelfand et al. [10]. To each matrix A g 
A^nxn(IHI), they associated functions, that we shall call quasicharacteristic 
functions, defined by 



f,j{X) = \Xld-A\,,, l<i,j<n, 



where | ■ \ij is the (z, j)— quasideterminant. Let us denote by A^'^ the submatrix 
of order (n — 1) obtained by suppressing the row i and the column j in the 
matrix A G A^„xn(H). Then quasideterminants are defined inductively by the 
formula 



\A\ij — aij 7 ^ ajqilA^'^ \pq) 



-1, 



J-PJ, 



where the sum is taken over the indices p,q & {!,...,«} with p ^ i,q ^ j , such 
that the quasideterminant of lower order jA^'-'Ipq is defined and it is non-null 
(see Proposition 1.5 of [5]). 
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When the matrix A is invertible, the entries of the inverse matrix are 
exactly a*^ = l^lj^"^- In the commutative case this gives the weh known for- 
mula a*^ = (—1)*+^ det j4'''Y det A. For quaternionic matrices, the norm of the 
quasideterminant \A\ij of ^ G A^„xn(IH) verifies 



Plyl •Sdet(^*'J) = Sdet(A). 
This is a particular case of Jacobi identity for quaternions (Proposition [14] 



(4) 



Remark 18. From Equation ([4]) it follows that the roots of the quasicharac- 
teristic functions are left eigenvalues. However none of those functions gives 
the complete spectrum, as shown in the next Example. Also notice that the 
definition of noncommutative left eigenvalue considered in |10[ subsection 8.2, 
p. 128] does not correspond to the notion we are discussing here. 



Example 19. Let A 



1 

k 

-3i 




j 

2k 



characteristic functions arc 



Then (7i{A) = {i,j,k}. The quasi- 



/ii(A) 
/12(A) 

/13(A) 
/22(A) 

/23(A) 

/32(A) 
/33(A) 

while /21(A) and /31(A) are not defined 



A-i, 

-(A-i)k(A-j), 

-(A-i) (3i-2k(A-j)-ik)"'(A-k), 
A-j, 

-i(A-j)k(A-k), 

-2k-(A-k)k(A-j), 
A-k, 



3.3 Characteristic map 

We now introduce the notion of a characteristic map whose roots are the left 
eigenvalues, thus generalizing the usual characteristic polynomial in the real and 
complex settings. As we shall see this notion fits naturally with the equation of 
order 2 given by Wood [55], as well as with the procedure proposed by So [53] 
in order to compute the left eigenvalues of 3 x 3 matrices. 

Definition 20. The map /j, : H H is a characteristic map of the matrix 
A S Mnxni^) if, up to a constant, its norm verifies ]/x(A)] — Sdet{A — Aid) for 
all A e H. 

Example 21. Let D = diag{qi , . . . , Qn) be a diagonal matrix. Then ^(A) = 
((71 — A) • • • {qn — A) is a characteristic map for D. Analogously for a triangular 
matrix. 
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Let us start with the 2x2 matrix A 



If A is a diagonal matrix, 



then ai{A) reduces to the elements in the diagonal. Otherwise, we can always 
suppose that 5 7^ (see Reniark[211). Moreover, Sdet(^ — Aid) does not change 
after elementary transformations (Corollary [TT]), for instance 



Sdet 



a — A 
c 



b 

d-X 



Sdet 







c - {d - X)b-\a - X) 



b 

d~X 



Then, as pointed out by Wood, computing the left spectrum is equivalent to 
finding the roots of a characteristic map like 



/i(A) =c-(d-A)6-i(a-A). 



(5) 



Remark 22. Huang [12 proposed another map when c 7^ 0, namely (A — 
a)c~^{X — d) — b. This polynomial is obtained by adding (A — a)c^^ by the second 
row to the first row. This expression is equivalent to b— {a — X)c~'^{d — X), which 
is the one given by Wood at the end of [3S] (there is a misprint in the original 
paper) . 

The left spectrum is not invariant by similarity. However, we shall use the 
following fact: 

Proposition 23. // P is an invertible real matrix then Sdet(y4 — Aid) = 
Sdet(P^P~^ — Aid). Hence A and PAP~^ have the same characteristic maps. 

Remark 24. Let A be a matrix of order n > 2, let Pais be the real matrix 
obtained by interchanging the rows a and /? in the identity matrix /„. Left (resp. 
right) multiplication by the matrix Pap switches two rows (resp. columns) of A. 
Now let i ^ j he two indices and let tt be any permutation of {1, . . . , n} sending 
i to 1 and j to n. Then tt can be written as a composition of transpositions, 
so by taking the product P of the corresponding matrices PajS we obtain the 
matrix PAP~^ where the initial entry aij of A has moved to the place (l,n). 



4 SPECTRUM OF MATRICES OF ORDER 2 

In the next paragraphs we shall classify the different possible spectra of 2 x 2 
quaternionic matrices depending on the rank of the differential fi^\ of a charac- 
teristic map. 



4.1 Preliminaries 

The characteristic map /i : H — ^ H given in ([5|) can be extended to a continuous 
(or even differentiable) map ^: S"' S''* on the sphere S"* = H U {oo}, because 
lim |^(A)| = oo when |A| — >■ oo. A rigorous proof of the following result can be 
found in Eilenberg-Steenrod's book 6, pp. 304-310]: 

Proposition 25. A polynomial map like fi and the power map X^ are homo- 
topic, hence they have the same topological degree, which equals 2. 
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From Lemma H] it follows that the differential of ^ is given by 
^l,x{X) = Xh-\a ~\) + {d- x)b-^x. 



(6) 



4.2 Classification of left spectra 

Now we are in a position to reformulate the following result from Huang and So 
[T5] . see also [HI [5]. Let A — ^ be a quaternionic matrix with & ^ 0, and 
denote 

ao = —b^^c, fli — b^^{a — d), A = — 4ao. 

Theorem 26 ([13]). The matrix A has one, two or infinite left eigenvalues. 
The latter case is equivalent to the following conditions: ao,ai are real numbers 
such that ao ^ and A < 0. 

Remark 27. We shall call spherical the infinite case, because the spectrum 
ai{A) = {(l/2)(a + d + bq): q"^ ^ A} is diffcomorphic to the sphere 5"^ C Hq = 
(ij,k). 

Let A S H be an eigenvalue of A, that is /i(A) = for the map /i in ([5]). 
In the next Propositions we shall apply Theorem |6] to the differential E = fi^\ 
computed in ([6]). Accordingly to the notation of Sylvester equation in Section 
12.31 we denote 

P={d-X)b-\ Q = b-\a-X). 
First we study the two non-generic cases. 

Proposition 28. If rarikfi^x ~ 0, then ao,ai are real numbers and A = 0. 
Moreover X equals (a + d)/2 and this is the only left eigenvalue of the matrix. 

Proof. We know from Theorem [6] that P = t G M. and Q = —t, then ai = —2t 
and 2X = a + d. From /i(A) = it follows that ao = +t'^, then A = 0. Now it 
is easy to check (using for instance Theorem 2.3 in [T3]) that A = a + <6 is the 
only left eigenvalue oi A. □ 

Lemma 29. Let A, B be two similar quaternions that do not commute. Then 
the equation A^ — (A + B)X + AB = has the unique solution X = B . 

Proof. If A ^ B is a solution, it follows from (A — i3)A = A{X — B) that A and A 
are similar, then 3?(A) — 3?(A) = 5ft(i3) and |A| = \A\ = \B\. By substituting in 
the equation we see that the real parts and norms dissapear, so we can suppose 
that A, _B, A are pure imaginary quaternions with norm 1. Hence A^ = — 1 = B^ 
so the equation reduces to {A + B)X = {A + B)B, which implies X = B, a 
contradiction. 

Alternatively, the uniqueness of A can be proved by using Theorem 2.3, case 
4.1, of the solution of quadratic equations in 14 . □ 

Proposition 30. //rank/i*A = 2 two things may happen: 
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1. either the spectrum is spherical and all the eigenvalues have rank 2; 



2. or the matrix has just one eigenvalue. 

Proof. By using the diffeomorphism a-\-hai{A') — iyi{A) we can substitute A by 

the so-called "companion matrix" A' = ^ ^ . Since the rank is 2, we 

[-flo -aij 

have from Theorem [S] that P = t + a and Q = —t + /3 with a,/3 S Hq — (i, j, k), 
\a\ = |/3| ^ 0. Then ai — —2t + /3 — a. The first possibility is that /3 = a, then 
ai = -2t. It follows from ^(A) = that ao = t^ + |/3|2 ^ and A = -4|/3|2 < 0. 
Then we have the spherical case. In particular A = (— ai — 2(3) /2. The other 
eigenvalues have the form (— oi +q)/2 with = — 4|/3|^, then the differential of 
II verifies P ~ t — q^^/2 and Q = —t — q/2, and so they have rank 2 too. The 
second possibility is that (3 ^ a. Then oi = —2t + f3 — a, oq = {t + a){t — f3), 
and Lemma shows that the only eigenvalue is A = t — /3. □ 

Now we consider the generic case. 

Proposition 31. //rankyU*A = 4 then the matrix has two different eigenvalues. 

Proof. Since the differential has maximal rank at the eigenvalue A, the matrix A 
cannot correspond to Propositions [28] or |30l hence all its eigenvalues are of rank 
4. Then by the inverse function theorem the fiber /x~^(0) is discrete (in fact 
compact) and its cardinal equals (Theorem [1]) the degree of the map /z, which 
is 2 by Proposition [55l Notice that the Jacobian is nonnegative by formula 

o. □ 

Remark 32. In [TB], Janovska and Opfer show that for quaternionic polyno- 
mials there are several types of zeros accordingly to the rank of some real 4x4 
matrix, but their procedure does not seem to have an immediate geometrical 
meaning. 



5 CHARACTERISTIC MAPS OF 3 x 3 MA- 
TRICES 

The only known results about the left spectrum of 3 x 3 matrices are due to 
So [24], who did a case by case study, depending on some relationships among 
the entries of the matrix. He showed that when n = 3 left eigenvalues could 
be found by solving quaternionic polynomials of degree not greater than 3. In 
general there is not any known method for solving the resulting equations. 

In the following paragraphs we shall develop a method for matrices of order 
3 which is analogous to that of Section l4] that is we shall find a map fj,A such 
that |/xa(A)| = Sdet(A — Aid). This time, however, the characteristic map ^ia 
will be in most cases a rational function instead of a polynomial (the latter 
occurs when the matrix A has some null entry outside the diagonal). 

a b c 

Let us consider the quaternionic matrix A — 



eX3x3(H). 
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5.1 Polynomial case 

We start studying the simplest situation, when there exists some zero entry 
outside the diagonal 

First, suppose that the matrix has the zero entry c — 0, that is 



Sdet(A-AId) = Sdet 



a - A b 
/ h 
p q r — X 



There are three possibihties: 

1. a b, h ~ 0, we have a triangular matrix, so we can take 
MA) = (r-A)(.g-A)(a-A); 



(7) 



2. if 6 = but h ^ 0, then Proposition [T2] allows us to reduce to the 2x2 
case and we obtain 



t^iX) = {q - (r ~ \)h~\g ~ \)) (a - A); 



(8) 



3. finally, if 6 ^ 0, we can proceed as follows. We create a zero in the first 
row by substracting to the first row Ci the multiple C2b~^{a — A) of the 
second column: 



Sdet(A- Aid) = Sdet 



b 

f ^ (g - X)b-Ha - X) g-X 
p ~ qb ^{a ~ X) q 




h 

r - A 



then we permute the two last columns in order to reduce to the 2x2 case. 
In this way we can take as a characteristic map the polynomial of degree 
3 



^iiX) =p- qb-\a - A) - (r - X)h-' (/ - {g - X)b-\a - A)) 



(9) 



Theorem 33. // the matrix A e A^3x3(H) has some zero entry outside the 
diagonal, then A admits a polynomial characteristic map. 

Proof. Let the entry be Oij = 0, with i j. Then accordingly to Remark [24] 
there is a real invertible matrix P such that the transformation PAP~^ does 
not change the characteristic maps and gives a matrix with 013 = 0. □ 

5.2 Rational case 

In the more general situation, when c ^ 0, we can compute the Study's deter- 
minant of the matrix A by creating zeroes in the first row. Then 



c 

/ — hc^^a g — hc^^b h 
p — rc^^a q — rc 



Sdet (A) = Sdet 

From Lemma [T^ and the results for 2 x 2 matrices it follows: 
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Proposition 34. If 0, then Sdet(^) is given: 

1. when g — hc^^b ^ 0, by 

\c\ ■ \g - hc-^b\ ■ \p - rc-^a - {q ~ rc-^b){g ~ hc-^byHf ~ hc-^a)\; 

2. otherwise, by 

\c\-\q-rc-'b\-\f~hc-'a\. 

Definition 35. We shall call pole of the matrix A G A^3x3(IHI) the point tta = 
g — hc^^b. 

Notice that tta is the quasideterminant |>l^'^|2i (see page [7]). 
By applying Proposition [34] to the matrix ^ — A Id we obtain the following 
characteristic map for A (we omit the term |c|). 

Proposition 36. Let A be a matrix of order 3x3 such that c 0. A charac- 
teristic map can be defined as follows: 

1. if TTA = g — hc^^b is the pole of A, then 

K'^a) = (g - (r - TrA)c^^b) (/ - hc~^{a ~ tta)) ; 

2. for A 7^ -KA we define 

fi(X) = {tta -X)[p-{r- X)c-\a - A) - (10) 
{q~{r- X)c~^b) {tta - A)-i (/ - hc-^{a - A))] . 

Remark 37. The map in (fTU)) is exactly the same formula given by So in [2^ 
p. 563], even if our method is completely different. This is why we chose to 
compute determinants starting from the right bottom corner. 

5.3 Continuity at the pole 



Up to now we have defined maps which verify |/i(A)| = Sdet(yl. — Aid) in norm. 
Since Sdet is a continuous map we have lim |/^(A)| = |/i(7r/i)|. However, the 

following example shows that /i may not be continuous at the pole tta- 

The pole tt^ = — i is not a 



i 1 

3i - k 1 

k -1+j + k 



Example 38. Let A = 

left eigenvalue; in fact 

h{tta) = (-1 + j + k + l)(3i - k - i) = 1 - i + 2j - 2k 
We observe that the directional limits 

limi4-i + eq) = -q{i + k)q-\2i-k) 

depend on g G H, hence lim /i(A) does not exist. 

X^TTA 
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Theorem 39. The characteristic rational map is continuous if and only if 
the pole TTA is a left eigenvalue of A. 

Proof. Assume that tta is a left eigenvalue. Let be a sequence converging 

to TTA- Then \^{qn)\ — Sdet{A — g„ Id) converges to Sdet(v4 — tta Id) — 0, that 
is /i(g„) = fJ.{TTA)- 

Now we shall prove the converse. The map /i defined in Proposition [33] is of 
the form 

H{X} = {nA-X)[p{X)-qiX)i7TA-X)-'f{X)], X^TTA, (11) 

while ij-{tta) — QiT^A)fiT^A)- Let as assume that lim /^(A) exists and equals 

niiTA)- We must check that /i(7r^) = 0. If /(tt^) = we have finished. Other- 
wise we deduce from pT]) that 

lim {X-TrA)q{X){X~7rAr' = -qiTTA). (12) 

From Lemma HDl it follows that the limit on the left side equals (/(tt^i), hence 
9(i'yi) = 0, which ends the proof. □ 

Lemma 40. Let Q = Q{X) be a continuous map and suppose that there exists 
the limit Iq — lim XQ{X)X^^ . Then Iq = (5(0). 

A — ^0 

Proof. Take a sequence of real numbers (e„)„ 0. Then 

^0 = lim£„(3(e„)e~^ = lim(3(e„) = (5(0). 

□ 

Notice that the differentiability of fj, at the pole tta is not ensured. 
It is an open question whether it is always possible or not to find a polyno- 
mial, or at least a continuous characteristic map for a given matrix A. 

5.4 Extension to the infinite point 

Each of the characteristic maps we have introduced up to now can be extended 
to the sphere 5"* H U {oo}. 

Proposition 41. The polynomial maps fi defined in Subsection \5. 1\ verify that 
lim |/x(A)| — oo. 

I A|— ^oo 

Proof. When 6 = 0, the result follows from formulae ([7]) and ([8]); when b ^ 0, 
the expression of is that of formula so 

Ia^(A)| ^ |(r - X)h-\g - X)b-\a - A)| \-p + gb-^a - A) + (r - X)h-^f\ 
|AP - |AP |A|2 

□ 
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Proposition 42. The rational map fi defined in the c^O case by formula 17^) 
can be extended to S'^ = MU {00} (maybe with a discontinuity at the pole tta)- 

Proof. We have 

/^(A) = {tta - A)p2(A) - (tta - A)gi(A)(7r^ - A)-Vi(A) 
with polynomials 

P2(A) - p-{r-\)c-\a-X}, (13) 
gi(A) = q-{r-X)c-'b, 
/i(A) = f~hc-'ia~X). 

Then 

KA)| > |(^A-A)p2(A)|-|gi(A)/i(A)| 

and 



MA)k |(^A-A)p2(A)| _ , ,,Ir-A| |a-A| 



1™ IMS ^ 1™ TTn = 1™ |c I |\| iM 1" r 

A|^oo |A|-^ |A|^oo lAI'^ |A|^oo |A| |A| 

□ 



5.5 Discontinuous case 

Let us now assume that the characteristic map 11 a defined in Section [36] is not 
continuous at the pole tta, or equivalently that tta is not a left eigenvalue of 
A (Theorem [39]) ■ Then the matrix i? = A — tt^ Id is invertible and its pole is 
ttb = 0. Moreover ai{A) = <Ji{B) +tia- On the other hand, from Proposition [T71 
we know that the spectra of B and B^^ are diffeomorphic, because ai{B^^) = 
ai{B)-\ 

Theorem 43. Let A be a matrix such that the pole tta is not a left eigenvalue. 
Let B = A — TT^Id. Then the matrix B~^ has a polynomial characteristic map. 

Proof. Accordingly to formula (|4]) the norm of the entry (1,3) of the matrix 
B~^ equals 

Sdet(B3'i)/Sdct(S) = IttbI/ Sdct(B) = 0, 
then Theorem [33l applies. □ 
Here is an alternative proof of Theorem Let m_b(A) — — Ai?(A), with 

i?(A)-p(A)+<z(A)A-V(A), A^O, 

be the characteristic map given in p6p (we assume ttb = 0.) Then it is imme- 
diate that A-R(A^^)A is a polynomial in A (of degree 3 and independent term 
—c~^) and we only have to apply the following result. 
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Proposition 44. Let ^ib he a characteristic map of the invertible matrix B = 
A — TTA Id, with a discontinuity at the pole ttb — 0. Then 

is a characteristic map for B^^. 
Proof. From Proposition 1101 we have 

Sdet(A"^ Id) • Sdct(B-i - Aid) • SdGt(B) = Sdet(A-i Id - B) 

then 

|A-3| • Sdct(S-i - Aid) • Sdet(S) ^ \hb{X'^)\- 

□ 

Remark 45. The idea that a rational map like R{X) can be converted into a 
polynomial A^^i?(A)A~^ with variable A~^ is due to So (see [24l Lemma 3.5, 
p. 558]). 

6 TOPOLOGICAL STUDY OF THE 3x3 CASE 

We shall consider separately the polynomial, rational continuous and discontin- 
uous cases considered in Section [5] 

6.1 Polynomial case 

Let us start with matrices having some null entry outside the diagonal (as we 
have seen this case can be reduced to the case c = 0). We know that the 
characteristic map is a polynomial of degree 3 which can be extended in a 
continuous way to the sphere S'^. Then since there is a unique term of higher 
degree 3, the map /i^ is homotopic to A'^, so it has topological degree 3. 

Proposition 46. Let X be a left eigenvalue of the matrix A with c = 0. Then 
the differential of the polynomial characteristic map ha in Subsection \ 5.1\ is 
given by 

1. if b,h = then 

/i*A(^) = -^(5 - A)(a - A) - (r - X){g - X)X - [r - X)X{a - A); 

2. ifb = 0,h^O then 

Xh-^g - A)(a - A) - (g - (r - X)h-\g - A)) X + (r - X)h-^X{a - 

3. otherwise, 

fi,x{X) = {qb-^ - {r - X)h-\g - X)b-') X 

+Xh-^ (/ - (g - X)b-^{a - A)) - (r - X)h-^Xb-\a - A). 
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The proof is a direct application of the derivation rules given in Lemma |31 
The expressions obtained are of the form PX + XQ + RXS = 0, whose rank 
can be computed with the method given in Section 12.41 

0" 

be a triangular matrix. The differential 



Example 47. Let A = 



a 

f 9 

[p q 

lJ,f:\{X) of the characteristic map 







/i(A) = (r-A)(.g-A)(a-A), 

at the eigenvalues \ — a, g,r is given, respectively, by (a — r){g — a)X, (g — 
r)X(a — g) and {r — g){a — r)X . Hence, unlike the case n = 2, the rank depends 
on the multiplicity of each eigenvalue, and can be either or 4. 



6.2 Rational case 

When none of the entries outside the diagonal is zero, the characteristic map is 
a rational function, with a distinguished point tta- 

Let us first suppose that the pole tta is a left eigenvalue. We know from (|10p 
that fiA is a continuous map on S"* of the form 

{nA-X)[p{X)-q{X){7rA-X)-'f{X)]. 

By examining formulae it is clear that p{X) is homotopic to — A^ by the 
homotopy 

tp - (tr - X){1 - t + tc^^){ta - A), < G [0, 1]. 

Analogously q{X) ^ X and /(A) ~ A, so ^a{X) is homotopic to 

{tta - X) [-A2 - X{^A - X)-^X] 

(notice that this map is continuous at A = 0), which in turn is homotopic to 
A'^ — A^ by the homotopy 

it^KA - A) [-A^ - X{tTTA - X)-^X] . 

Finally the homotopy A^(A — t) shows that the map ^a is homotopic to A'^. All 
these homotopies can be extended to the infinity. 
Hence we have proved 

Proposition 48. When the rational characteristic map fiA 'is continuous it has 
topological degree 3. 

On the other hand, suppose that tta is not a left eigenvalue. Then the 
polynomial case applies to (A — tta Id)~^ by Theorem H51 So we do not have to 
use the local theory of degree, whose main difficulty is the need of considering 
homotopies which are admissible with respect to the domain 51 of definition, see 
D p. 28]. 

Corollary 49. Any 3x3 quaternionic matrix has at least one left eigenvalue. 

Proof. In all cases the eigenvalues (or its inverses) can be computed as the roots 
of a continuous map /i of degree 3; then, /i~^(0) is not void (see Section[5]). □ 
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6.3 Final remarks 



In order to simplify the computation of the rank, by taking B = A — tt^ Id we 
can always assume that the pole is ng — 0. 

Proposition 50. For a 3 x 3 matrix with c ^ the differential of the charac- 
teristic map II given in formula \1(J\) at the point A ^ ttb = is 

lJ-*xiX) = 

X [-P +{r- A)c-i(a -\) + {q-{r- \)c-^h){-X)-^{f - hc'^a - A))] 
+ {-X)Xc-^{a - A) - {-X)Xc-^b{-X)-\f - hc'^{a - A)) 
-{-X){q - (r - X)c-^b){-X)-^X{-X)-\f - hc-\a - A)) 
+ [(-A)(r - A)c-i - {-X){q -{r- X)c-^h){-\)-^hc-^] X. 

When is continuous but not differentiable at the pole tt^ the rank at 
IT A could be computed by taking a different characteristic map fipAP-^- In 
particular, by moving around the entries outside the diagonal of a given matrix 
(see Remark [23] of Section one can obtain up to six different characteristic 
maps, each one with a different pole. 

It is an open question whether there exists a matrix A verifying that all the 
poles of the matrices PAP~^, with P real, are eigenvalues. If such an example 
does not exist then the non-polynomial case would not be necessary. 



7 EXAMPLES 

Here are some miscellaneous examples. 

Example 51. Discontinuous map. Let A be the matrix given in Example [38l 
Then 

i i 1 

B A - TT^ Id 3i - k i 1 

k -1+j+k i_ 

is an invertible matrix with pole tt^ = 0. By computing the quasideterminants 
we obtain the inverse matrix 

4i-2k -4i + 2k 

-l-3i + 8j-6k l + 3i-3j + k -5j - 5k . 
ll + i-8j-8k -l-i + 3j + 3k -5j + 5k 



B-^^1- 
10 



Its polynomial characteristic map is 



(14) 



10 - Ai - 2iA - ^iA(2i - k)A - i^A(l + j + 2k)A - ^A(j + k)A(2i - k)A. 

On the other hand, the rational characteristic map for B is 

^b(A) = -Ai?(A) = -A [1 + k + lA + Ai - A^ + (j + k + Ai)A~i(2i - k + A)] 
and one can check that Ai?(A^^)A equals (up to a constant) the map ([T4| . 
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Example 52. An eigenvalue of rank 3. Let 



A 



j 1 

2i -k 1 

2 i 2j I j + k -i k 



The characteristic map is 

= 2 - i - 2j + (1 + j - k) (j - A) + (i + k + A) (21 + (k + A)(j - A)) , 
For the left eigenvalue A = the differential is 

^,o(^) = kX + Xi + (i + k)Xj, 



whose real associated matrix M 



"0 


-2 














-2 

















2 





-2 






has rank 3. 



Example 53. A matrix which can not be reduced to the polynomial case c = 0. 
We shall exhibit a matrix A G 7^3x3(111) such that for any real invertible matrix 
P all entries in the matrix PAP^^ outside the diagonal are not null. Let 
A = T + iX+jy + kZ, with T,X,Y,Z £ X3x3(R). Then 



PAP-^ = PTP-^ 



iPXP-^ +jPYP-^ + kPZP- 



which means that the matrix PAP^^ has a null entry if and only if the same 
happens for the real matrices PTP^^, PXP~^, PYP~^ and PZP~^. Moreover, 
from Remark [24l we can suppose that the null entry is at the place (1, 3). Now, 
recall that A is the matrix associated to a linear map — > H'^ with respect 
to the canonical basis 61,62,63 and that the real matrix P ~ (Pij) represents 
the change to another basis vi,V2,V3, that is Vj — J2iP^''^i where P^^ = (p'-'). 
Assume that T, X are the matrices associated to two rotations R"^ M"^ with 
the same rotation axis C and rotation angles +7r/2 and — 7r/2 respectively. Then 
the nullity of the entry (1,3) means that Tv^, Xv^ G (w2, W3), which implies that 
either ^3 is in the direction of the axis, i.e. C = (W3), or it is orthogonal to the 
axis, in which case (u2, W3) = C'^ (otherwise it is impossible that v^, Tv^ and Xv^ 
lie in the same plane) . Now, suppose that Y, Z are two other rotations with axis 
C and rotation angles ±7r/2. If C and C are different and not perpendicular, 
then it is impossible that Yv^,Zv^ G {v2,v^)- 

Example 54. A continuous rational characteristic map. The pole tt^ = 1 + j 
is an eigenvalue of the matrix 

"0 -j i' 

p q r 
with p, g, r G H arbitrary, p,q =/= 0. In fact, 

Mtt^) = (q - (r - 1 - j)k) (-1 + j + j(-l - j)) = 0. 
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Example 55. Generic polynomial case. Let 



■ k 0' 
3i-j -i i 
1 - 2k j -j 

The characteristic map is /Lt(A) = (— 1— k+Ai)A(k — A) hence <Ji{A) = {k, 0, — i — 
j}. The differential of fi at each eigenvalue is 

M*k(^) = (-l-i + k)X, 
M*oW = -(l + k)Xk, 
= X(l + 2i + k). 

Then the matrix A has three different eigenvalues, all of them with maximal 
rank. 

Example 56. Two eigenvalues, one of null rank, the other one of maximal 
rank. Let 

F-i-j 
k i i 

1-i j 

This time ii{X) = (1 + k - Ai)A(i + j + A) so ai{A) ^ {0, -i - j}. We have 

M*oW = (l + k)X(i+j), 
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